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We use the semiclassical approach combined with the scaling results for the diffusion coefficient to 
consider the two-level correlation function J?(e) for a disordered electron system in the crossover 
region, characterized by the appearance of a macroscopic correlation or localization length, ^, that 
diverges at the metal-insulator transition. We show new critical statistics, characterized by a non- 
trivial asymptotic behavior of J?(e), to emerge on both sides of the transition at higher energies, 
and to expand to all energies larger than mean level spacing when ^ exceeds the system size. 
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Single-electron spectra in disordered metals are gov- 
erned [1-3] by the universal Wigner-Dyson (WD) statis- 
tics [4] which are applicable to a large variety of different 
quantum systems. The principal conjecture is that spec- 
tra of complicated systems are statistically equivalent to 
the eigenvalues of the random matrix hamiltonians [5] 
constrained due to the presence of some general symme- 
tries. The main feature of the WD statistics is the level 
repulsion at all energy scales. For a disordered system, 
it is due to the electronic eigenstates being extended. 

With the increase of disorder, the system undergoes 
the Anderson metal-insulator transition [6]. At the other 
side of the transition, in the insulating phase, there is 
no correlation between the energy levels of the localized 
states, and the spectral statistics prove to be Poisson. 

As it has become clear recently, there should exist the 
third universal level statistics, applicable exactly at the 
transition point. Its existence has been first asserted by 
Shklovskii et al. [7] who have suggested that the nearest- 
level probability distribution, P(e), was a universal hy- 
brid of the WD one at small e and the Poisson one at 
large e. Systematic analytical studies of the universal 
critical statistics have been started in Ref. [8]. Namely, 
the two-level correlation function, has been studied in 
the framework of the diagrammatic approach and the 
scaling hypothesis, and found to be drastically different 
from both known universal statistics. In contrast to the 
Poisson statistics, the long-range level correlations are 
still present, but they are different from the WD ones. 
The mapping of the analytical results of Ref. [8] onto the 
one-dimensional plasma model with the power-law inter- 
action [9] illustrated the presence of the level repulsion 
at large energies. It shows itself in P(e) having the large 
energy asymptotic behavior intermediate between Gaus- 
sian for the WD statistics and exponential for the Poisson 
ones [9], rather than being just Poisson as in Ref. [7]. 

All the three universal statistics are exactly applica- 
ble in the thermodynamic limit, Lq, oo. The critical 



statistics are applicable within an energy band with a 
fixed (but arbitrary large) number of the levels, centered 
exactly at the mobility edge, = £c, while the WD and 
Poisson statistics are applicable for energy bands with 
£g > £c (the metallic region) and < £c (the insulating 
region), respectively. However, when the sample size Lq 
is finite, there is a smooth crossover between the metal- 
lic and the insulating phase. In this Letter, we present 
the analytical description of an appropriate crossover be- 
tween the WD and Poisson statistics. 

The crossover region is characterized by the appear- 
ance, on approaching the Anderson transition, of a new 
macroscopic scale, ^, which is the correlation or local- 
ization length at the metallic or insulating side of the 
transition, respectively. It diverges at the mobility edge 
as [6] 

^ = i^\l-gjgcr, (1) 

where is some microscopic length (that can be of the 
order of the elastic scattering length, i = v^t), is the 
conductance at the scale , gc is its critical value, and 
i> is the critical exponent that depends on the dimen- 
sionality, d, and the universality class. We will show the 
new critical statistics to emerge at large energies in the 
crossover region, when <^ < Lq, and to become uni- 
versal and applicable to all energies [10] e ^ A (with A 
being the mean level spacing) near the transition, when 

Calculations in the crossover region may be performed 
rigorously within the diagrammatic approach, similar to 
that used at the mobility edge [8]. In this Letter we 
will use instead a simplified description based on a more 
intuitive and illuminating way of treating the spectral 
correlations developed recently by Argaman, Imry, and 
Smilansky [11] within the semiclassical approach. 

We consider the two-level correlation function (TLCF): 
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Here p{E) is the electron density of states at the energy 
E for a particular realization of disorder, (...) stand for 
averaging over all the realizations, p = {p[E)), and A = 
1/[pLq). The spectral form factor is defined as 
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A transparent semiclassical expression for it, based on 
the Gutzwiller trace formula [12], has been obtained by 
Berry [13]. Combining the Berry's expression with a gen- 
eralized Hannay and Ozorio de Almeida's sum rule [14], 
Argaman, Imry, and Smilansky [11] have shown that 
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where /3 = 1, 2, or 4 is for the Dyson unitary, orthogonal, 
and symplectic ensembles, respectively, p{t) is the av- 
eraged classical probability density to perform periodic 
motion of period i at a given energy [15]. It is differ- 
ent from the averaged classical return probability for a 
diffusing electron (which is contributed to by any path 
returning the point of origin, including those with mis- 
matching initial and final momenta) simply by a constant 
factor. 

One identifies two opposite time scales: ergodic, re- 
quired for the diffusion motion to fill the entire phase 
space, t ^ Terg = Lq/D, and diffusive, t ^ Terg, with D 
being the diffusion coefficient. (We involve here neither 
ballistic, t ^ t, nor quantum, i > h/A, regimes). They 
correspond to the energy scales s <^ g and s > respec- 
tively, where g = Ec/A = hpDL^'^ is the conductance 
in units ofe^ /%, and Ec = h/Tg^g is the Thouless energy. 

Efetov has shown [2] the statistics in the ergodic regime 
to be the same as in the Wigner-Dyson random matrix 
theory [5]. There is no dependence either on d, or on 
D, as only the spatially homogeneous diffusion mode (i.e. 
the q = Fourier component of the diffusion propagator) 
contributes to the correlation function. It corresponds to 
the saturation of the probability density [13,11] in Eq. 
(4), p(i) = const. Then K(t)cxt results in i?erg(s) ~ 1/*'^, 
and thus in a logarithmic spectral rigidity [16]. 

At the shorter times, t ^ Terg, the probability distri- 
bution of a diffusing electron is a standard Gaussian: 
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The return probability p(i) is obtained by setting r = r^. 
Substituting it into Eq. (4), one obtains in the diffusive 
regime [11,15]: 
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Performing the Fourier transform, one reproduces the di- 
agrammatic result of Altshuler and Shklovskii [3]: 
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where is a positive constant for d > 2, C3 ~ 1. Al- 
though -Rdif(s) is not universal, it depends only on the 
conductance g and the dimensionality d so that the dif- 
fusive regime is governed by the one-parameter scaling. 
Note that it corresponds to the levels attraction [17]. 

We can still use such a semiclassical approach to de- 
scribe the spectral statistics in the crossover region and 
even at the mobility edge. The trick is to substitute 
the diffusion constant in Eq. (6) by the scale-dependent 
diffusion coefficient known from the scaling theory of lo- 
calization [6]. 

In the crossover region, <^ C ^ Lq, the natural en- 
ergy scale related to ^, Eq. (1), is 
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that is the mean level spacing within a correlation vol- 
ume. It is known from the scaling theory of localiza- 
tion [6] that at scale L ^ ^ the conductance shows an 
Ohmic behavior, g = aL'^~'^, with a = gci^'"^. The elec- 
tron propagation at this scale remains diffusive, with the 
scale-independent diffusion constant 
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The corresponding time scale is given by i > ^ jD(^ ~ 
h/gcAp^ ~ h/A^, as ~ 1 for d = 3. At this time scale, 
the TLCF is still given by Eq. (7) as the rescaling of D 
is absorbed by an appropriate rescaling of g. 

Totally new statistics emerge at the larger energies [10], 
in the critical regime, 



s > A|/A 



t < h/A^ , 
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corresponding to the diffusion at the scale A ^ ^. In this 
regime, the conductance is almost scale-independent and 
the diffusion is anomalous, with the coefficient given by 
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To the first approximation, ffcr(A) = gc, and the time 
dependence of the anomalous diffusion may be found by 
combining Eq. (11) and A'^ = Z)cr(A)i which leads to the 
well-known [18] result: 



i3cr(i)~^0(i) = (ffc/V)'^' + 
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Substituting this to Eq. (6), we find the form factor to 
be time-independent that results in 



Ro{s)=0, s>A|/A. 
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Such an absence of the correlations in the critical regime 
of Eq. (10) results from using the self-consistent approx- 
imation (12). For A ~ io ^, it corresponds to the 
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absence of correlations at the mobility edge found in the 
same approximation diagrammatically [8]. (To prove the 
diagrammatic cancellation, it was necessary to use not 
only the scaling relations (12) but also the causality re- 
quirements which are quite subtle in the diagrammatic 
technique but automatically taken into account in the 
present formalism). It means that the self-consistent ap- 
proach is insufficient to describe the spectral statistics in 
the critical regime. One must find Z)cr(A), Eq. (11), more 
accurately, allowing for the scaling dependence of gcr(A). 

The standard scaling equation for the conductance, lin- 
earized near the critical point g = Qci has the form [6] 
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dlnA 



gc 



(14) 



Integrating this on the metallic side with the initial con- 
dition g[A = ig) = g^, one obtains in the critical regime 



ffcr(A) = ffc 1 + (A/0 
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with ^ defined by Eq. (1). It yields, after substituting into 
Eq. (11), an additional time-dependence of the diffusion 
coefficient, on top of that in Eq. (12), 
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Dcr{t) = Do{t) [l + (ffciA|/a)- 
thus leading to the time- dependent form factor. 
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where a is a numerical coefficient of order 1. This time- 
dependence leads to the following result for the TLCF in 
the critical regime of Eq. (10): 
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with Adjs being a numerical coefficient depending only on 
the dimensionality and the universality class. Therefore, 
there are three regimes in the crossover region: ergodic 
for t > Terg (including the quantum times, t > h/A) 
where statistics are Wigner-Dyson, diffusive for Terg ^ 
t > h/A^ where statistics are the same as in the non- 
ergodic regime in metals, Eq. (7), and critical, Eq. (10), 
where the new statistics emerge [19] described by the 
TLCF (18). When ^ increases in approaching the tran- 
sition, Eq. (1), the critical region (10) is expanding (see 
Fig. 1). Finally, at ^ ~ Lq, both A| and Ec become of 
the order of A, so that the WD regime shrinks to the 
quantum limit (s <C 1 O i ^ h/ A), the diffusive regime 
disappears entirely, and the critical regime expands to 
the whole [10] region w > A (Fig. 1). There the TLCF 
becomes completely universal: 
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FIG. 1. Schematic phase diagram for the statistics in the 
crossover region from the metall (where the correlation length 
^ ^ ^o) to the mobility edge (the dotted line, ^ ^ Lo) to 
the insulator (where the localization length ^ l^). Rescaled 
WD regime is reduced to the WD regime at ^ ^ Lo and to the 
Poisson regime at HLo 0, see Eq. (21). Here St = H/tA, 
and = Ep/A. 

This expression is not changed with a further increase of 
^, as for ^ > Lq diffusion in the whole sample is anoma- 
lous and qualitatively the same as for ^ ~ Lq- Equation 
(19) coincides with that obtained directly at the mobility 
edge [8] within the diagrammatic approach. 

In the above procedure, it seemed to be straightfor- 
ward to estimate the numerical coefficient a in Eq. (17), 
and thus Adp in Eq. (18). However, we could find neither 
its value, nor even its sign, as this procedure was rather 
oversimplified. There are two reasons for this. First, the 
semi-classical expression (4) has been shown in Ref. [11] 
to be equivalent to the two-diffuson (and two-cooperon) 
diagrams of Ref. [3]. However, at the mobility edge and 
thus also in the crossover region, there are infinitely many 
of relevant diagrams. All of them have been proved to 
give parametrically the same contributions at the mobil- 
ity edge [8]. Then, to get a parametrically correct answer 
it is sufficient to consider just the simplest typical contri- 
bution that can be described in the semiclassical language 
of Eq. (4). Second, we used the Gaussian form, Eq. (5), 
of the probability distribution P[r,t), albeit with the 
time-dependent diffusion coefficient D[t). In fact, in the 
region of anomalous diffusion, an exact expression for P 
is not known. It remains causal, however. This, together 
with the scaling representation of its Fourier transform at 
the mobility edge, P(g,e) oc [F[q^, A^/Qq''' — ie]~^ (with 
Ae ~ [pe)~^l'^ being a typical diffusion length during 
the time h/e and F being an arbitrary sc&Ymg function), 
allowed to prove [8] that the TLCF, Eq. (19), remains 
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unchanged, except for the numerical coefficient. Such 
a proof that involves all possible diagrams is extended 
straightforwardly to the crossover region, as will be pub- 
lished elsewhere. 

Finally, let us describe the statistics in the crossover 
region, ^ ^ Lq, on the insulating side of the tran- 

sition, where ^ is the localization length. The level cor- 
relations are only due to the states confined to the same 
localization volume. The TLCF is obtained then as a su- 
perposition of independent contributions, R^, from each 
volume 

i?insW^tfW + ^i?|(s^) , (20) 

where the 6[s) function describes the self-correlation of 
the levels. 

At smaller energies, e < A|, the argument of is 
smaller than 1, so that it should be qualitatively the same 
as the Wigner-Dyson TLCF which, in turn, becomes al- 
most a constant for s <C 1. Thus, we obtain in the region 

i?i„3(^) ^ S{s) + A R^^(^s « S{s) - A . (21) 

The regular part is small as only the small fraction of 
the levels, (^/io)"^ = A/A|, are correlated. The Poisson 
statistics emerge in the thermodynamic limit, where this 
regular tail vanishes. 

At larger energies, we arrive at the critical regime of 
Eq. (10), where the scaling relations differ from those on 
the metallic side, Eqs. (15) and (16), only by the sign 
of the ^ dependent corrections. Thus, the TLCF on the 
insulating side may differ from that on the metallic side, 
Eq. (18), only by a numerical coefficient. Although this 
coefficient is unknown, as explained above, it should be 
equal to that in Eq. (19). Indeed, there is no distinction 
between metal and insulator at the critical point where 
A = A|, and Eq. (20) should coincide for s ^ 1 with 
Eq. (19), obtained at the mobility edge. Then one should 
substitute -Rme for R^ into Eq. (20), which then becomes 
identical to Eq. (18), obtained in the critical regime (10) 
on the metallic side. So, the critical statistics, universal 
at the mobility edge, are equally applicable on both sides 
of the metal-insulator transition in the crossover region. 

ACKNOWLEDGMENTS 

We are thankful to B. L. Altshuler, Y. Gefen, V. V. 
Lebedev, A. D. Mirlin, and B. Shapiro for usefull discus- 
sions. A.G.A. is grateful to GIF for financial support and 
ICTP for kind hospitality. V.E.K. and I.V.L. gratefully 
acknowledge travelling support under the EEC grant No. 
SSC-CT90-0020. 



[1] L. P. Gor'kov and G. M. Eliashberg, Zh. Eksp. Teor. Fiz. 

48, 1407 (1965) [Sov. Phys. JETP 21, 940 (1965)]. 
[2] K. B. Efetov, Adv. Phys. 32, 53 (1983). 
[3] B. L. Altshuler and B. I. Shklovskii, Zh. Eksp. Teor. Fiz. 

91, 220 (1986) [Sov. Phys. JETP 64, 127 (1986)]. 
[4] E. P. Wigner, Proc. Cambridge Philos. Soc. 47, 790 

(1951); F. J. Dyson, J. Math. Phys. 3, 140 (1962). 
[5] M. L. Mehta, Random matrices (Academic Press, Boston, 

1991). 

[6] P. A. Lee and T. V. Ramakrishnan, Rev. Mod. Phys. 57, 
287 (1985). 

[7] B. I. Shklovskii, B. Shapiro, B. R. Sears, P. Lambrianides, 

and H. B. Shore, Phys. Rev. B 47, 11487 (1993). 
[8] V. E. Kravtsov, I. V. Lerner, B. L. Altshuler, and A. G. 

Aronov, Phys. Rev. Lett. 72, 888 (1994). 
[9] A. G. Aronov, V. E. Kravtsov, and I. V. Lerner, Pis'ma 
v ZhETF 59, 39 (1994) [JETP Letters 59, 40 (1994)]. 

[10] We do not consider here the ballistic energy scale, 
e > ft/r ~ Ep. Statistics of the electrons at the appro- 
priate time scale, t ^ t, should be insensitive to the 
Anderson transition and remain in the crossover region 
basically the same as those in the metallic phase stud- 
ied by Altland and Gefen [20]. Although this regime is 
of no interest for the level-level correlations in the criti- 
cal regime, as the TLCF vanishes in the thermodynamic 
limit, it still could be important for integral statistical 
characteristics, such as the variance of the levels number 
in a given energy band, and might contribute to the sum 
rule discussed in Ref. [8]. 

[11] N. Argaman, Y. Imry, and U. Smilansky, Phys. Rev. B 
47, 4440 (1993). 

[12] M. C. Gutzwiller, J. Math. Phys. 12, 343 (1971). 

[13] M. V. Berry, Proc. R. Soc. Ser. A 400, 229 (1985). 

[14] J. A. Hannay and A. M. Ozorio de Almeida, J. Phys. A 
17, 3429 (1984). 

[15] We consider the range of energies where the correlation 
function (2) and thus the form factor (4) do not depend 
on E so that Eq. (4) has a less general form than in Ref. 
[11]. The coefficient of proportionality in Eq. (4) differs 
from that in Ref. [11] because of the additional factor, 
l//0^, in our definition of R{s), Eq. (2). 

[16] To obtain not only the envelope but the correct oscil- 
latory behavior of R{s), more subtle treatment of the 
longer, quantum times, t ^ ft/A is required [13]. Here we 
always consider energy scale e ^ A that corresponds to 
t < ft/A. 

[17] R. A. Jalabert, J.-L. Pichard, and C. W. J. Beenakker, 

Europhys. Lett 24, 1 (1993). 
[18] B. Shapiro and E. Abrahams, Phys. Rev. B 24, 4889 

(1981); Y. Imry, Y. Gefen, and D. J. Bergman, Phys. 

Rev. B 26, 3436 (1982). 
[19] Equation (7) and (18) are parametrically matching at 

the boundary of the diffusive and critical regimes, i.e. at 

s ~ Af/A, both yielding l/s. Indeed, it results from Eq. 

(8) that g oc {A^/Af-^''^ for L > ^ > 4, thus giving 

the prefactor (Af/A)^-^''^ in Eq. (7). 
[20] A. Altland and Y. Gefen, Phys. Rev. Lett. 71, 3339 

(1993). 



4 



